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Abstract. This work describes one parametrization of quantum states and channels 
and several of its possible applications. This parametrization works in any dimen- 
| sion and there is an explicit algorithm which produces it. Included in the list of 

■ applications are a simple characterization of pure states, an explicit formula for one 

£NJ | additive entropic quantity which does not require knowledge of eigenvalues, and an 

(— i ■ algorithm which finds one Kraus operator representation for a quantum operation 

^ ' without recourse to eigenvalue and eigenvector calculations. 
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> \ 1. Introduction 

t> 

The interest in quantum information processing has brought added attention to 
questions of parametrizations of positive matrices. Indeed, at least two foundational 
ingredients in the theory of quantum information, viz., quantum states and quantum 
channels involve positive matrices. See, for instance, Therefore it appears to 

be of interest to obtain descriptions (parametrizations) of the set of positive matrices. 

The purpose of this note is to provide such an explicit parametrization by adapting 
similar results on positive definite kernels in This parametrization is rather intri- 
cate (and nonlinear), but in this paper we show some benefits of its use: the parameters 
of pure states are easily described; the parameters of a tensor product can be deduced 
from those of its factors; at least one entropic quantity (to be described in Section 3.3) 
is easy to compute; the Peres- Horodecki criterion for separability can be described ex- 
plicitly described as inequalities for these parameters in low dimensions; purifications 
of qubits can be explicitly parametrized, and finally one Kraus operator representa- 
tion for quantum operations can be computed without requiring any knowledge of the 
eigenvalues/eigenvectors of the associated Choi matrix. 

The key features of this parametrization which are worth emphasizing are that it 
works in any dimension (indeed, it applies to matrices whose entries are operators 
themselves); there are explicit formulae for the parameters in terms of the entries of 
the matrix; while these are rather intricate, there is one computationally attractive 
algorithm which produces this parametrization ; several quantities of interest can be 
computed via these parameters; the algorithm exploits the fact that the inherent struc- 
ture in the matrix is inherited by its Schur complements and this results in the algorithm 
yielding (at least) one Cholesky decomposition of the matrix. Indeed, this last feature 
is precisely the reason why a Kraus operator representation of a quantum channel is 
produced by these parameters, without any need for eigenvalues or eigenvectors. 
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This paper is organized as follows. After preliminaries concerning definitions, the 
main result on the parametrization of positive semidefinite matrices is stated and il- 
lustrated via some examples in the next section. We deal with states in Section 3 and 
with quantum channels in Section 4. An Appendix contains the sketch of a proof of 
Theorem 1 on the parametrization of positive matrices and describes how to view a 
positive matrix as a matrix with a so-called displacement structure and then describes 
an algorithm for the computation of the parameters. 

2. Preliminaries 

In this section we introduce terminology and state the main result concerning the 
parametrization of positive (semidefinite) matrices. 

2.1. Quantum states and channels. The state of a rf-dimensional quantum system is 
described by a d x d positive density matrix of trace 1, that is, a positive element of 
trace 1 in the algebra Aid of complex d x d matrices. States described by rank one 
density matrices are called pure states. 

A quantum channel is a completely positive map $ : A — » C(7i) from a C7*-algebra 
A into the set £(H) of all bounded linear operators on the Hilbert space Tt (in the 
situations most frequently met in quantum information processing, A = Aid and 
C(TC) = Ai d ', while $ is also required to be trace preserving). By the Stinespring 
theorem, [THj, Theorem 4.1, such a map is the compression of a *-homomorphism. For 
A = Aid, there is a somewhat more explicit representation, given in |2| ( see also 0)- 
Thus, $ : Aid — > C(7i) is completely positive if and only if the matrix 

(2-1) S = S*= ME k ,)} d k , J=l 

is positive, where E k j, k,j = l,...,d, are the standard matrix units of Aid- Each E k j 
is a d x d matrix consisting of 1 in the (k,j)th entry and zeros elsewhere. We notice 
that if X = [Xk,j] k j =1 , then Y = [Ykj] k j =l = $(X) is given by the relations 

(2.2) Y k j = <&(Ei jm ) k jXi }Jn . 

l,m 

This shows that there is a one-to-one correspondence between the set of completely 
positive maps onAlj with values in C(7i) and the set of positive matrices in Ai d®C-(H). 

For a linear map $ : Aid — ► Aid the adjoint <f> is defined with respect to the Hilbert 
space structure on Aid given by the Hilbert-Schmidt inner product (A, B) = Tr(A*B), 
where A* denotes the adjoint of A. It is easily seen that $ is trace preserving if and 
only if l> is unital ($(/) = I). 

2.2. A parametrization of positive matrices. We describe a parametrization of the 
positive matrices in Aid ® £(H), with TC allowed to be infinite-dimensional. Note that 
if dim(H) = 1, such matrices are precisely dxd positive matrices with complex entries. 
To that end, some elements of dilation theory ([HI) are needed. Let Hi and H2 be two 
Hilbert spaces, not necessarily finite dimensional, and let C(Hi,TC2) denote the set of 
all bounded linear maps operators from 7ii into 7i 2 . The operator T e £(7ii,7i 2 ) is 
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called a contraction if ||T|| < 1. The defect operator of T is Dt = (I — T*T) 1 / 2 , where 
T* denotes the adjoint operator, (as well as the complex conjugate in case T is just a 
complex number) and let T>t denote the closure of the range of Dt- To any contraction 
T E C(Hi, H 2 ) is associated the unitary operator U (T) : Hi © V T * — * ri 2 © V T by the 
formula: 

(2.3) U(T) = 

Now let {J- k } k= i be a family of Hilbert spaces and consider {T k j \ k, j = 1, . . . , d, k < j} 
a family of contractions such that (i) T k;k — for k — 1, . . . , d and (ii) for k < j, 

(2-4) T k , j eC{V Tk+l .,V Tl ._ i ). 

Unitary operators U k j are associated to this family by the recursions: U k ,k — the 
identity operator on jF fc , k = 1, . . . , d, while for j > k, 

Uk,j = Uj-k(Xk,k+i)Uj-k(X k)k+2 ) . . . Uj_k(Xk,j)(Uk+ij © iv v * )) 

where Uj^k(Xk,k+i) denotes the unitary operator defined from the space 
(®£i2>r* +M+ J © (Vr k+hk+l ®V ri k+l ) © (& m=l+1 V rk k+m ) 
onto the space 

(©^i^r fc+1 , fc+m ) © {Vri^ ®V Tk>k+l ) © {®l =l+1 V rk k+m ) 
by the formula 

^_ fc (r fc)fc+ ,) = / ®u(r kfk+l ) © /. 

The foregoing considerations from dilation theory will now be applied to matrices. 
Let S = [Sk,j] k j =1 be a matrix such that S k ,j € £(H) and let S k , k — L* k k Lk± be a 
factorization of S ktk , k = l,...,d. Denote the closure of the range of L kjk by Tk-, 
k = 1, . . . , d. For a family of contractions T k j as in the previous paragraph, denote by 
R k j the row contraction 

R k ,j = [ r* jfc+ i, D r * ktk+ r kM2 , D rik+i . . . D r%j _T k j ] 

and by C kj j the column contraction 

c k,j = Tj-2,jD rj _ ld , T kij D Tk+1 . . . . D Vj _ ltj ] , 

where "f stands for matrix transpose. We now obtain the following characterization 
and structure of a positive matrix. 

Theorem 1. The matrix S = [S k j] k j =1 as above, satisfying S* k = S k j, is positive if and 
only if i) S k k > 0, k — 1, . . . , d and ii) there exists a family {Tkj \ k, j — 1, . . . , d, k < j} 
of contractions such that T ktk = for k — 1, . . . , d and Equation \2-4\j is valid, and 

(2.5) S kJ = L* ktk (R kij -iU k+ i tj -iC k+l j + D r * k+i . . . I)y ; A'k.jlh-, . , • • • Ih . , )/-.,,,• 
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For a proof see Appendix 1. The Tkj,j 7^ k form the parametrization of positive 
matrices referred to in Section 1. Note the are just some fake parameters, included 
in the statement of the theorem to avoid an artificial separation of the j — k + 1 
case from that for other values of j Note that when the dimension of TC = 1, these 
contractions, Tkj, are complex numbers in the closed unit disc) . 

Though this result looks rather intricate it provides a true parametrization of the 
set of positive matrices, has a useful physical interpretation, and does not depend on 
the dimension of 7i. Furthermore, by adopting certain natural conventions, this can 
be even turned into a one-one parametrization of quantum states. One illustration of 
this is provided in Section 3, for the d = 2 case. 

Let us next look at some concrete examples of Theorem 1. For d = 2 the result 



is well-known (see, e.g. 



For d = 3, the structure is more interesting, even in 

be a positive matrix. 







' S n 


Sl2 


<Sl3 


the scalar case dimH = 1. 


Thus, let S = 


Q* 
D 12 


S22 


5*23 






°13 


D 23 


S33 


Theorem ^ gives: 


S12 = L^T 


L2-^22, 







23 



-^22^23-^33, 



s 



13 



L{i (ri 2 r 2 3 + -Dr* 2 ri3-Dr 23 ) £33, 



and we can notice the connection between the structure of S and spherical geometry 
in the formula for S13. Thus, for Ti 2 = cos6*, T 2 3 = cos#i and Ti 3 = cos0, the above 
formula for S13 reduces to the cosine law in spherical geometry. 

We also write the case d = 4 explicitly, with 7i of arbitrary dimension; in this case, 



if S = [S, 



k,j\k,j=V 



then 



5, 



k,k+2 



L 



k.k 



Sk,k+1 — L* kk Tk,k+lLk+l,k+l 

r Vfe+if '^+1^+2 + Drt . Ti 



1,2,3, 



L 



k+2,k+2i 



1,2, 



and 



S 
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-^n (ri2r23r34 + -Dr* 2 ri3-Dr 23 r34 + ri2-Dr* 3 r24-Dr 34 
--Dq 2 ri3r2 3 r 2 4.Dr34 + d v *^d v * 3 Ti±d v ^d v ^ l 44 . 

Lattice Structures and Time-Dependent String Models It is worthwhile to digress briefly 
to point out that in the Appendix, positive matrices are shown to be a special case 
of matrices with displacement structure, and that this results in an algorithm for the 
calculation of the parameters r^. Furthermore, due to the special choice of a lower 
triangular F(—t) (in the definition of a displacement structure) made in the appen- 
dix, this algorithm has the so-called lattice-structure. The same structure produces a 
Cholesky factorization of S, 

S = (®k=l L k,k) G l,d Gl , d {®k=l L k,k) , 
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where the upper-triangular operator G^d is defined recursively by Gk,k 
1, . . . , d, and for 1 < k < j < d, 



I, k 



(2.6) 



G 



k,0 







Uk,j-\Ckj 
D Th ...D T .. 



These lattice structures can alternatively be described in the form of a time-varying, 
discrete transmission-line (string) For more details see j3], where this connection is 
described in the context of displacement structures. In the invariant time case, this 
transmission line is familiar in marine seismology, where the one-sided perfect reflector 
is given by the interface air- water, see jSj- The implications of this classical model 
for quantum states remain to be fully worked out. However, this model explains the 
intriguing fact that the number of summands in the formula for Sk,k+h / = 0, 1, . . ., in 
Theorem 1, is precisely the Ith Catalan number. These observations will be explained 
elsewhere. 



3. Parametrization of states 

In this section we parametrize (finite dimensional) quantum states by using The- 
orem ^ and we show several applications of this parametrization. Throughout this 
section, except in Subsection 3.2, we will write the Tkj and the Dv* kl of the previous 
section as g k ,i and dk,i respectively, so as to distinguish states from channels. The 
explicit formulae for the parameters, gk,i, depend on the basis of Aid and it is often 
convenient to use selfadjoint bases containing the identity (see, for instance, [12] )• 
Such an orthogonal basis can be obtained from {Ekj}ij=i as follows: 

ft = E k>j + E j>k , k < j, 

ft = - (E jtk - E kJ ) , k > j, 
i 

h d 1 = I d , h{ = /»*- 1 ©0, I <k < d, hj=J-^—^ (ht^il-d)). 
For d = 2, we deduce 

hi = h, fu = ci> fix = °~2, hj = o- 3 , 

where a\, 02, 03 are the Pauli matrices. In the Pauli basis a hermitian matrix with 
trace 1 will be representated as p = + Y^k=x fikO~k) and its positivity is equivalent 
to the Bloch sphere condition 1 — ||/3|| 2 > 0. By using Theorem d the positivity of p 
is equivalent to a cylinder condition on the parameter (@3,g): 

Iftl < 1, 

Px-iP2 = (l-\f3 3 \ 2 ) 1/2 g, \g\<l, 

and in order to ensure a one-to-one parametrization we choose g = if j /?3 1 = 1. 
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For d = 3 we obtain the Gell-Mann matrices: 





1 








' 3 " V3 


1 








hf = I 3 , hl = 





-1 








1 






















-2 



A hermitian matrix of trace 1 will be represented as 

P 



By using Theorem ^ the positivity of p is equivalent to the conditions: 
l + & + ^|>0, l-/3 2 + A> , l-^>0, 

712 - *721 = ((1 + ^|) 2 - /5 2 2 ) 1/ V2, |<7l 2 | < 1, 

v l — /5 2 H 7=) 7 (1 7=) ' 923, 1 023 1 < 1, 



723 — ^732 



V3' 



713-^731 



(l-/?2+-^) 1/2 (l-^|) 1/2 (gvtito + (1 - I^i2| 2 ) 1/2 (1 - Itol 2 ) 1 ^) , |<?i 3 | < 1. 



The case d = 4 was explicitly written inj^j, but with respect to the basis a k ® a r In 
any dimension we can write 



P 



V 1=2 k^j / 



and the parameters f3i appearing on the diagonal of p satisfy a system of linear inequal- 
ities, while the parameters gkj associated to the off-diagonal elements by Theorem ^ 
are almost independent: they belong to the closed unit disk and satisfy the boundary 
conditions (|2.4j) . 

3.1. Pure states. It was already noticed in [S] that the purity of qubits can be explicitily 
checked via the parameters of Theorem ^ Here we extend this result to any dimension. 
Let p be a d- dimensional state, 



P 



V 1=2 k^j J k,j 



■J 



We notice that 



and 



Pk,k = 1 + ^f3i(hf)k,k, k = l, 



,d, 



Pkj = lh,i ~ ijj,k, k < j. 

Let {gkj} be the parameters associated to the positive matrix p = [pk,j\'j. J=1 by Theo- 
rem ^ 
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Theorem 2. A state p is pure if and only if the parameters gkj are zero except for 
those indices for which Pk,kPj,j 7^ 0, in which case \gk,j\ = 1- 

Proof. Assume p is pure. This implies that each matrix 

Pk,k Pk,k+l 
_ P*k,k+1 Pk+l,k+l 

has rank at most 1. This happens in case either at least one of pk,k, Pk+i,k+i is zero or 
if |<7fc,fc+i| = 1. It remains to see what happens if pk t k 7^ 0, pij = for k < I < m and 
Pm,m 7^ 0. In this situation we must have p p j = for k < p < I < m, therefore g v ^\ = 

1/2 

for k < p < I < m, and so, pk, n , ~ " ' ' 



Pkk9k,mPm% by (|2.5|1 . Then, the fact that 



Pk,k Pk,m 
Pk,m Pm,m 



has rank at most 1 implies \gk,m\ — 1- 

Conversely, we deduce from (|2.5|) that p 
as follows: let i\ < . . . < i/. be the indices of the nonzero diagonal elements of p. Then 



-v v* where the vector v is described 



J PJ*2 



1/2 
Pi\,i\ i 

1/2 * 
Pii,i-$i\iii ' 



\ u PJik 

while the other entries of v p are zero. 



1/2 * 
Pih,ik9ii,i2 



□ 



Thus, the vector in C d representing the pure state can also be written down in terms 
of the parameters. It was noticed in |2j that these parameters can be used in producing 
a one-one parametrization of purifications of qubits. 

3.2. Tensor products. We next discuss producing parameters for tensor products of 
positive matrices, in terms of the parameters of the factors entering the product. 

1 A 12 



Let us begin with a very simple example. Take S\ 



A 



12 



and S2 



B 



12 



B\2 

1 



Then 



1 


-Bl2 


A 12 


A 12 B 


r>* 
B X2 


1 


A 12 B^ 2 


A 12 


A* 


By2.A\ 2 


1 


B\ 2 


R* A* 

- D 12 /1 12 


A* 




1 



s 1 ® s 2 



(note that we use the so-called right Kronecker product representation of the tensor 
product, see pQ). Formula ()2.5|) gives: 



r 12 — Bi 2 , r 2 3 — Ai 2 B^ 2 , r 34 — b 12 , 
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^12(l-|£l2| 2 ) 1/2 

(l-|A 12 | 2 |5 12 p)i/2' 

A l2 {\ - \B 12 \^ 
(l-|A 12 | 2 |i? 12 | 2 )V2' 

The fact that Ti4 has a simpler structure than expected indicates that there is 
additional structure that can be explored fruitfully. It turns out that in this case it is 
more convenient to search for a block lattice structure first. 

Let S x = [Slj] d kJ=1 , S 2 = [Sjy fc J=1 be two positive matrices. Also let {T^}, {r 2 ^} 
be the corresponding parameters associated by Theorem^ to Si, respectively, S* 2 . For 
a matrix T we use the notation T® p = T ® . . . ® T , to denote its p-fold direct sum. 

p times 

Theorem 3. The positive matrix Si ® S 2 is a d x d block matrix with structure given 
by (EH) with: 




Lk,k — G\ d ,(S\ k )® d , k — 1, . . . , d, 

r fc j = (ry ed ', i<k<j<d, 

where G\ d , is the upper triangular Cholesky factor of S* 2 . 
Proof. First, it is easily checked that 

s^s^iGi^rksijrt (Gi. d ,r 

and then, some algebra with formula ()2.5|) gives that the parameters associated by 
Theorem [T] to the positive matrix [(-S'fcj) ed ] fc -_ x are F^- = {l}j)® M ■ □ 

It is quite convenient to use the Cholesky factor of S 2 since it can be explicitly 
computed in terms of the parameters T 2 ,^, as shown by (|2.6j) . 

3.3. Entropy. In this section we consider an entropy-like number E(p) that can be 
explicitly computed in terms of the parameters {gkj} of a state p. This quantity, 
E{p) = ^Tr (logp), is mentioned in [T7j, as one of a list of candidates for a quantum 
notion of entropy. Further material about quantum entropies is summarized in [17J. 

One simple way to motivate E(p) is to start with the classical Kulback-Leibler in- 
formation number, 

D(f\\g) = [ f(x)log(f(x)/g(x))dx, 



where / and g are probability densities and x = {x\, . . . , x n ). If / and g are Gaussian 
with covariance matrices P, respectively, R, then 

D(f\\9) = -~ (logdetP/T 1 + tr(J - PR- 1 )) . 



By setting, for instance i? = ^, we see that —E(p) is upto a constant D(f\\g). Thus, 
loosely speaking, in this view the classical analogue of a state is a Gaussian whose 
covariance matrix the density matrix is. The foregoing suggests the consideration of 
the following entropy of a d- dimensional state, 



k=l 

where {Xk}f =1 is the set of eigenvalues of p. This is a well-known formula giving 
a functional of entropy type (see [H]), and also plays a significant role in convex 
optimization, [T3j . 

E behaves quite differently from the von Neumann entropy, for instance, pure states 
have vanishing von Neumann entropy while E(p) = — oo for a pure state p (however, 
see below for one variation on E(p) which remedies this). Still, E has some interesting 
properties, two of them described by the following result. 

Theorem 4. (a) Let p be a state with parameters {gkj}- Then 

1 / d 

(3-1) E(p) = - £>gp fc , fc + £>g(l - \g hj \< 

\fc=l k<j 

(6) Let p and ip be two states. Then 

E(p®^) = E(p)+E(iP). 

Proof, (a) Using Algorithm 10 and Th 11 yields (cf., Th 1.5.10 in 

det P = (n^. fc )n( 1- w 2 )- 
\fc=i y k<j 

Taking the logarithm in this formula we obtain 1)3.1)) . 

(6) This follows easily from the formula for the determinant of a tensor product of 
matrices. □ 

We also notice that E(p) < logo?, with equality for p n = ... = p dd = i which 
implies the following maximum entropy principle: if p is a given state, then for any 
dimension d' there is a state ipo of dimension d! such that 

max E(p ®ip) — E(ip ), 

i> 

where the maximum is taken over all states ip of dimension d 1 . 
One small variation on E is provided by the formula: 

1 ' 
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where A^ , j = 1, . . . , I, are the nonzero eigenvalues of p. We see that Eq is still additive 
and p is pure if and only if E (p) = 0. However, to apply 1)3.1)1 . one must first find the 
restriction of p to its support. 

The key utility of using the parameters g^j is that there is no need for any eigenvalue 
computations for finding E(p). Contrary to the situation with eigenvalues, the g k j can 
be related to the entries of p via explicit formulae. 

3.4- The Peres- Horodecki Criterion. The parameters {gk,j} of a state p can be also 
used to explicitly write the finite set of inequalities characterizing the separability of 
2x2 and 2x3 states. For 2x2 states we have the following result. 

Theorem 5. Let p be a 2 x 2 state with parameters {gkj}- Then p is separable if and 
only if: 

1/2 1/2 1/2 1/2 , 777^ 
P22 P33 + Pll P44 «12«13«24«34 ^ 

Pll P44 1012023034 + ^125 , 13^23fl'34 + 012^23024^34 ~ ^12013023024^34 1 

and the system of inequalities 

dl2 {d 2 3 + (1 - l^l 4 l) 1/2 ) > 1012023 -9vihu\, 

d 3 4 (c?23 + (1 - |^l 4 l) 1/2 ) > 1034023 ~ 034^14|, 

1/2 1/2 . 1/2 1/2 , ,,2 lxl/2/i 17 2 |\1/2j \ 

P22 P33 +Pll P44 ^12(1 - \h 13 \) 1 (1 - |/i 24 |) 1 d u > 

p\lp4A 1^12^14^34 + ^12^13(1 - \hl 4 \) 1/2 g u + 0i 2 (l - \hl 4 \) 1/2 h 2 4d U - ^12^13^14^24^341 

admits solutions h\±, his, ^24 *w i/ie closed unit disk, subject to the boundary condition 
(J23D, where d kJ = (1 - |^| 2 ) 1/2 . 

Proof. The result is a consequence of the Peres-Horodecki characterization of 2 x 2 
separable states, jEllIH], and Theorem [T] □ 

Similar inequalities can be written for 2x3 separable states. 

4. Quantum channels 

In this section we analyse some consequences of Theorem for the structure of 
quantum channels. We can immediately exemplify the case of binary channels. 

Example 6. A detailed analysis of quantum binary channels is given in [22 • We show 
here how Theorem relates to that analysis. It is showed in ^31 that any quantum 
binary channel $ has a representation 

*(A) = U[$ tA (VAV*)]U*, 

where U,V G U (2) and $ t ,A has the matrix representation 



T = 



1 











h 


Ai 








*2 





A 2 





h 








A3 
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with respect to the Pauli basis {/, 01, a 2 , c 3 } of M. 2 . See |2Z| also for similar normal 
forms. We can obtain (formula (26) in [22]) that 



S, 



$t.A 



l+£ 3 + A 3 tx-it 2 Ai + A 2 

£i + it 2 l — t 3 — A 3 Ai - A 2 

Ai - A 2 1 + £ 3 - A 3 ti - it 2 

Ai + A 2 ti+it 2 l-£ 3 + A 3 



Similarly, by formula (27) in [2*2*] . 



5, 



3>t,A 



1 + tg + A 3 ti + it 2 Ai + A 2 

1 + £ 3 - A 3 Ai - A 2 h + it 2 

ti - it 2 Ai - A 2 1 - £ 3 - A 3 

Ai + A 2 h + it 2 1 - £ 3 + A 3 



It is slightly more convenient to deal with S = [Sk,j]ij = i = 2S|, t . Formula (|2.5|) gives: 

= 1 + £ 3 + A 3 ; >S 22 = 1 + £ 3 — A 3 ; 

S33 — 1 ~~ ^3 — A 3 ; S44 = 1 — £ 3 + A 3 ; 
Ti 2 = 0, T 34 = 0; 

'23 



C C 1 / 2 F oV 2 

093 — o 22 1 23 o 33 , 



so that 



Ai — Aj 



23 



(l + £3-A 3 ) 1 /2(l-£3-A 3 )V2^ 



513 — 'S'ii 2 ri3-Dr23'S'33 ■ 
(£ 1 + z£ 2 )(l+£ 3 -A 3 ) 1 / 2 



,1/2 



so that 



13 



((1 + £ 3 - A 3 )(l - £ 3 - A 3 ) - (Ai - A 2 ) 2 )V2(i + h + A 3 )V2 ' 



c c 1 / 2 n f c 1 / 2 

D 24 — 22 -t^rio 1 24>->44 ) 



so that 



Finally, 



(£ 1 +z£ 2 )(l-£ 3 -A 3 ) 1 / 2 



24 



((1 + £ 3 - A 3 )(l - £ 3 - A 3 ) - (Ax - A 2 ) 2 )V 2 (l - £ 3 + A 3 )V 2 • 



5*14 — S 



1/2, 
11 1 



■ri3r 23 r 24 + D r * XuDr 24 )sl{ 2 . 



For brevity, we omit writing out the formula for r 14 . We deduce that $ t ,A is completely 
positive if and only if the following eight inequalities hold: 

Sk,k>0, fc = l,...,4, 

|r 23 | < 1, [rial < 1, |r 24 | < 1, |r 14 | < 1. 

Further, we know what happens in the degenerate cases (i.e., the cases where any of 
these inequalities become equalities). Thus, the implication of Skk = for some k on 
the structure of $ t ,A is clear. Also, if |r 23 | = 1, then necessarily t\ = t 2 = and 
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A x + A 2 = (l + t 3 + A3) 1/2 ri4(l-t3 + A 3 ) 1/2 for some contraction r M . If either |r 13 1 = 1 
or | r 2 4 1 = 1, then necessarily T 14 = and = Sl( 2 (— T 1^2^24) S^ 2 . 

We notice that this result is of about the same nature as that in [22] • This is because 
the first step of (j2.5jl . specialized to 4 x 4 matrices ( viz., the j = k+1 step), is precisely 
Lemma 6 in [22] which is used for the analysis in [22] ■ If we had instead used the block 
version of ()2.5j) . i.e., viewing S*$ as a 2 x 2 matrix with entries themselves 2x2 matrices, 
then we would deduce precisely Theorem 1 of 22\. What we basically have done here 
is that we used ()2.5j) in order to deduce in a systematic way the condition that R$ t A 
in Theorem 1 of [22] is a contraction. One advantage of doing this is that it works in 
higher dimensions. 

It is noted that the correspondence between S$ and the parameters r&j is nonlinear. 
Only for the first step is the correspondence affine and therefore can be used in the 
analysis of extreme points in the case d — 2, as it was done in [22]- This seems to be 
unclear for d > 2, at this moment. □ 

4-1. Capacity. It is worth remarking that the suggested notion of quantum entropy 
from Section (3.3) could be used to posit a notion of channel capacity for quantum 
channels, which is trivially additive. Several numbers have been suggested to define 
the capacity of a quantum channel, with additivity as a desirable requirement. See, for 
instance, [231 12H]- While additivity of some these notions has been demonstrated for 
special classes of channels, such as entanglement breaking channels, conjectures in the 
direction of additivity for some other notions were recently disproved (|25j). Therefore, 
it might be useful to have some other possible candidates. Motivated by the discussion 
in Section 3.3 we introduce for a quantum channel $ : Ada —* M-d 1 the number 

£>($) = -I log detS*, 

where N = dd! . Trivially, D(§) is additive, and there is an explicit expression for it in 
terms of the parameters Y^j of S<$>: 

1 / N 

(4.i) £>(<&) = -- nriogs M +5> g (i - ir fcJ r 

\fc=l k<j 

(b) 

Once again, -D($) > logn and a minimum capacity principle holds: given $, there 
exists a quantum channel \I> such that 

minD($ ® 

is attained. As an example of explicit computation of D (<$>), we consider binary quan- 
tum channels as in Example [HI For simplicity, assume t\ = £2 = ^3 = 0. Then, using 
the parameters calculated in Example El we deduce 

£>($) = -1/4 (2 log +2 log 



+ log(l-|^! 2 )+log(l-|^| 2 )) 



12 



4-2. Connections with Kraus representations. In this subsection, we provide one ex- 
plicit Kraus operator representation for any channel, which does not require either 
eigenvalues or eigenvectors of the associated Choi matrix. While Kraus representa- 
tions are not unique, they are very useful in explicit computations involving quantum 
channels. For instance, one can write down at least one Stinespring dilation in terms 
of a Kraus representation, and thus at least one mock unitary operation for a quan- 
tum channel; a formula for entanglement fidelity of a channel can be computed in 
terms of them; sufficient conditions for either unitarity or the entanglement breaking 
property can be checked; they play a role in quantum error correcting codes, quantum 
tomography etc., See, for instance, [Tot UU\ I7j. 

Consider a quantum channel $ : Aid — ► Aid- A familiar representation of quantum 
channels is the Kraus representation, [14| . 

d? 



k=l 



where A k e Aid, k 



d 2 , are called the generators of $, and Ylt=i AkA k = I 



The connection between the Kraus representation and is given by the formula 
(4.2) 
where 



A 



S$ = A* A, 
row(Ai) 
. row(A d 2] 



and 

row(A fc ) = [ (A k ) n 



{A k ) ld {A k ) 2d 



k 2d 



(A 



k dl 



(A k ) dd ] ■ 



This relation shows that while the generators of a Kraus representation are not unique, 
they will give the same channel provided that (|4.2|) holds. In particular, we can choose 
Kraus generators using the Cholesky factorization of S$. As mentioned before, the 
algorithms for finding the parameters, T k j, systematically yield a Cholesky factorization 
of a positive matrix. We give the details for binary channels. 

Example 7. Let S$ : M.2 M-2 be a binary channel and let be the parameters 

of $ given by Theorem [TJ Using formula ()2.6|) . we deduce 



G 



with 
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X 







ZL 33 







Dr 12 L22 


(D ri2 T 23 - T* 12 T 13 D r23 


) -^33 XL 44 








Dr 13 D r23 L 33 


YL44 











A"l4 -^24^34-^44 


Dr 12 


T23T34 + 


A^a-^f 23^24-^34 





— rj 2 ri3_Dr2 3 r34 + r^riar^I^A^ — Tl 2 D r * 3 Ti4Dr 24 Dr 34 , 

13 



and 



Y - D Tl3 Dr 2 . i T 3 4 - Dr 13 T* 23 T 2 4Dr 34 - T* 13 T 14 Dr 24 D rs4 , 

z = r 12 r 2 3 + D r * 2 Ti 3 D r2:j , 
w = r 12 r23r 3 4 + D r » 2 r 13 D r23 r34 + r 12 -Dr* 3 r24-Dr34 
-D ri r 13 r* 3 r 24 Dr 3i + d vi d vi y u d v ^d Tza . 



The Kraus generators are 

Al = \ L > 

ZL 33 w i 





ri2-C'22 

WL 44 



A, 







-Dri 3 £^23-^33 ^-^44 



(^r 12 r 23 



.4, 





ri 2 ri3-Dr 2 3^ l 



Dr V2 L22 



33 



12 

XL 



11 



D Ti4 D T24 Dy 34 L 44 



The advantage is that, once the parameters of S$ associated by Theorem ^ are 
known, the Cholesky factorization can be computed by (|2.fij) . Indeed, the algorithms in 
the appendix for calculating the Tkj, systematically compute a Cholesky factorization. 
Also, the orthogonality properties of the Cholesky factorization might be useful in some 
situations. These constructions make sense in infinite dimensions, except that further 
qualifications are needed for the trace-preserving property to be meaningful. 



5. Appendix 

In this appendix, a sketch of the proof of Theorem is first given. This proof does not 
immediately yield a constructive procedure for finding the r^ . For that purpose, we 
next discuss briefly the notion of displacement structure (which works for a family of 
matrices with operator entries). A key feature of a displacement structure is that all 
possible Schur complements of such a family of matrices inherit a related displacement 
structure (this is what enables the determination of Cholesky factorizations of each 
member of the family). This is stated in Algorithm |H| We next show how any positive 
matrix can be imbedded into a family of matrices with displacement structure. Typi- 
cally, a family with displacement structure admits more than one such representation. 
We choose one such representation - the specific choice for F(—t) enables the associa- 
tion of a lattice structure to the corresponding version of Algorithm |H1 Now, if Ji is of 
finite dimension, this specific choice of displacement structure made for positive matri- 
ces, results in a significant simplification of Algorithm |S| This is stated in Algorithm 
In particular, a series of explicitly determined finite-dimensional contractions, 7fc(t), 
is produced by the algorithm. In terms of these, 7fc(t), there is a simple and explicit 
formula for the parameters, Tkj- Thus, the determination of the parameters, T^, is 
fully constructive when Ti is finite dimensional (much of the procedure survives even 
when Ti is infinite dimensional). 

5.1. Parametrization of positive matrices. We now sketch a proof of Theorem ^ Since 
the diagonal entries of S intervene in the parametrization only through the Lkk {Skk = 
L kk *L kk ), we assume, without loss of generality, that = for all k — 1, . . . , d. By 
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a result of Kolmogorov (Theorem 1.3.1 in [3]), there exist a Hilbert space IC containing 
Hand isometries V k G £(TC, IC), k = 1, . . . , d, such that 

(i) The set {V k TC \ k = 1, . . . d} is total in IC; 

(it) S k j = V k *Vj, k,j = l,...,d. 

Since V\ is an isometry, we can identify Ti with the range of V\ and assume V\ = P^/Ti., 
where P^ denotes the orthogonal projection of IC onto Ti. Since V 2 is an isometry, there 
is a Hilbert space V 2 and an operator D 2 G £(^2, such that 

Wi = [ V 2 D 2 ]:H®V 2 ^IC 

is a unitary operator and V2 = W\/Tl. Set K\ — K,, K 2 — H. ® T> 2 . By induction we 
obtain a family of unitary operators W\ G C(fC 2 , /Ci), . . ., W^-i G C(JCd, K,d-i) such 
that 

V fc = W 1 ...W fc _ 1 /W, k = 2,...,d. 

In particular, 

(5.1) S kd = P%W k ...W j - 1 /H, k< 3 . 

Now, we notice that by a suitable identification of spaces, 





Tl2 


a 


b 


Wi = 


Dr 12 




A 12 







A21 


A 22 



Multiplying this operator on the left by the unitary operator R(r* 2 ) © I, we obtain a 
unitary operator with matrix representation 



/ 


X 


Y 





A'n 


A' 





A' 


A' 



It follows that necessarily, X = and Y — 0. Iterating this construction and using 
(|5.1j) and (ii), we deduce (|2.5|) . 

5.^. Algorithms As discussed above, the notion of displacement structure is useful in 
producing an algorithmic procedure for finding the T^j- The systematic study of the 
displacement structure was initiated in ^U] . A main theme of the theory is the recursive 
factorization of matrices with implicit structure encoded by a so-called displacement 
equation of the form 

(5.2) R(t) - F(t)R(t-l)F(t)* = G(t)J(t)G(t)*, teZ. 

Given R(t — 1), we need to know {F(t),G(t), J(t)} in order to determine R(t) using 
(|5.2|) . We shall not seek to determine R(t) by explicitly applying (|5.2|) . Instead, we 
shall use the fact that R(t) is a "low rank" modification of R(t — 1) and exploit it to 
determine R(t) more efficiently This idea will operate as follows. Use as input data 
{F(t), G(t), J(t)} and the Cholesky factor of R(t — 1), say L(t — 1), and then compute 
the Cholesky factor of R(t) without determining R(t), 

(5.3) R(t) = L(t)L*(t) . 
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The following algorithm, see [Hj, tells us how to compute the columns of L(t) from 
the columns of L(t — 1) and knowledge of {F(t),G(t), J(t)}. Here we describe the 
algorithm for special types of {F(t),G(t), J(t)}, which subsume the case of positive- 
definite matrices R(t). 

Let lk(t) denote the nonzero part of the k—th column of Lit). Let also 1/a/ d k (t) 
denote the top entry of l k (t) and define 

hit) = y/drft)h(t). 

That is, the top entry of h(t) is normalized to 1, and we also obtain the equivalent 
triangular factorization for R{t), 



(5.4) 



R(t) = Lit) D- 1 it) L* it) 



where the diagonal entries of D{t) are the {dkit)} and the columns of L(t) are the 
{hit)}- Here, L(t) is a unit diagonal lower triangular matrix. Let R k {t) denote the 
Schur complement of R(t) with respect to its leading k x k block. Suppose Fit) is 
lower triangular with diagonal entries f k (t) and let F k {t) denote the matrix obtained 
by deleting the first k rows and columns of F{t). 

Algorithm 8. Assume we know 

{L(t-l),Dit-l),G(t),F(t),Jit)}, 

and that R{t) satisfies the displacement equation A5.ty) . Then, for each k, the Schur 
complements, Rkit), of R{t) with respect to its leading k x k block, satisfy similar 
displacement equations, 



(5.5) 



Rkit) - F k (t)R k (t - l)F*(t) = G k {t)J{t)G* k {t), 



where the G k {t), and the triangular factorization \5.Jj\) of R(t), can be obtained from 
the following recursive construction: 



(5.6) 



hit) 





Gk+i(t) 



[ F k (t)l k (t-1) G k (t) ] 



r k {t) h* k it)Q k (t) 
J(t)9t(t) Pt(t)®k(t) 



where g k (t) is the top row ofG k (i), and h k (t) and p k (i) are chosen so as to satisfy the 
relation 



fkit) g k (t) 
h k (t) Pk(t) 



dkit - 1) 
Jit) 



fkit) g k (t) 
h(t) p k (t) 





'd^t) 




Jit) 



It is shown in jH] that it is always possible to find h k (t) and p k {t) as above. Choices 
that result in array (lattice) form descriptions are also possible and are described in 
detail in [B]. 
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Displacement Structure for Positive Matrices Let S be positive. For t = 0, . . . , d 

I St+l,t+2 St+l,d 

S*+l,t+2 I $t+2,d 



R(-t) = 

c* 

L °t+l,a 

so that, R(0) = S$ and R(—d + 1' 



t times 




I. Then, for t = 0, 



d time 



F(-t) 





/ 



G(-t) 



and 







J{-t) = J 



t+l,t+2 



f+l,t+2 



t+l,d 





c* 













We easily check that 

(5.7) R(-t) - FR(-t-l)F* = G(-t)JG(-t)*, t=l,...,d-2, 

so that, d — 2 steps of Algorithm |H] will exploit the whole information encoded in 
R(Q) = S$. 

We notice that Algorithm|Hlis greatly simplified with the special choice of R(t), F(t), G(t), 
when Tt is of finite dimension. Thus, write guif) = [uk(t) Vk(t)] with respect to 
the decomposition of J and since fk(t) = 0, we deduce dk(t) = gk(t)Jgk(t)* = 
v>k(t)y>k(t)* — Vk(t)vk{t)* . The positivity of S$ is therefore tested by the condition 
v>k(t)uk(t)* — Vk{t)vk(t)* > for all t and all k. If this holds for some t and some k, 
then there exist explicitly determined contractions, 7&(£), satisfying such that Vk{t) = 
Wfc(i)7fc(t). In this situation we introduce 



I 

-ikity 



-ik{t) 
i 



(i- lk (t) lk (tyy 1/2 
o 



o 

(/-7*(*)*7*(*)) 



-1/2 



and then the generator recursion in Algorithm |H1 reduces to a simpler form. 

Algorithm 9. The generator recursion for the displacement equation (|5.7J) associated 
to a positive matrix S has the form: 



(5i 





Gk+i{t) 



FG k (t-i)e k (t 



I 




G k (t)e k (t) 




/ 



Algorithm IH1 produces an explicit connection between the parameters 7&(t) and the 
entries S k j of S. Namely, a direct calculation yields the following result. 
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Theorem 10. T kJ = 7j_ fc (l - k)* , 1 < k < j < d. 
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